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ABSTRACT 

We have vastly increased the C IV statistics at intermediate redshift by surveying the thousands 
of quasars in the Sloan Digital Sky Survey Data-Release 7. We visually verified over 16,000 C IV 
systems with 1.46 < z < 4.55 — a sample size that renders Poisson error negligible. Detailed Monte 
Carlo simulations show we are approximately 50% complete down to rest equivalent widths W T ~ 0.6 A. 
We analyzed the sample as a whole and in ten small redshift bins with approximately 1500 doublets 
each. The equivalent width frequency distributions f(W T ) were well modeled by an exponential, with 
little evolution in shape. In contrast with previous studies that modeled the frequency distribution as 
a single power law, the fitted exponential gives a finite mass density for the C IV ions. The co-moving 
line density dAfciv/dX evolved smoothly with redshift, increasing by a factor of 2.37 ± 0.09 from 
z = 4.55 to 1.96, then plateauing at dJVcw/dX « 0.34 for z — 1.96 to 1.46. Comparing our SDSS 
sample with z < 1 (ultraviolet) and z > 5 (infrared) surveys, we see an approximately 10-fold increase 
in dAfciv/dX over z « 6 -J 0, for W T > 0.6 A. This suggests a monotonic and significant increase in 
the enrichment of gas outside galaxies over the 12 Gyr lifetime of the universe. 

Subject headings: galaxies: halos - intergalactic medium - quasars: absorption lines - techniques: 
spectroscopic 

Online-only material: color figures, machine-readable tables 



1. INTRODUCTION 

The study of the large-scale structure of the universe 
provides top-level constraints on models of galaxy evolu- 
tion. Heavy elements are produced in the stars of galax- 
ies. A variety of feedback processes move these met- 
als from the sites of production into the intergalactic 
medium (IGM), enriching the material for future gener- 
ations of stars. The cosmic enrichment cycle generically 
refers to the movement of gas from inside galaxies to the 
IGM and, possibly, back again (perhaps many times). 
The amount of heavy elements, the number of ionizing 
photons from galaxies and quasars, and the spatial dis- 
tribution of material are driven by hierarchical structure 
formation and galactic processes such as star formation 
and feedback. Spectroscopic surveys of quasars yield a 
random sample of intervening absorbing gas clouds that 
can be used to constrain the on-going and summativc 
enrichment processes in the universe. 

C IV AA1548, 1550 doublets are important tr acers of 
the I GM and its evolution from z ~ 6 to today (|Steidel 
Barlow fc Tytler|19981JEllison et al |2000[|Songaila 



et all 120091 IBecker et al.||2009l ICooksey et alJ|2T3To] 
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D'Odorico et al. 2010 Simcoe et al. 
sition has been well-studied at 1.5 ^ z 
following reasons. First, it is a strong transition of a 
common metal. Second, it is observable outside the Lya 
forest, where it becomes easier to identify. Third, it red- 
shifts into optical passbands at z — 1.5. Lastly, it is a 
resonant doublet, which gives it distinctive characteris- 
tics that enable surveys to be largely automated. 

Observations of C IV doublets constrain the cumula- 
tive effect of the cosmic enrichment cycle. More specif- 
ically, the number and strength of C IV absorbers are 
affected by: the amount of carbon produced by all pre- 
vious generations of stars; the spatial distribution of the 
element, driven by feedback processes; and the total ion- 
izing radiation available to maintain the triply-ionized 
transition. 

The intermediate-redshift C IV studies have tradi- 
tionally found that the doublets follow a power law in 
the column density distribution function, wit h a slope 
of a w —1.8, throughout the redshift range ( Songaila 
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12 < logVciv < 15 absorbers ((Songaila) 2001 
|berg et al. 2003 Songaila||2005[) . Improve' 



_Boksen-[ 
observations, 

pushing to bluer wavelengths (and hence lower redshift), 
have shown that f^ciVj the C IV mass density relative 
to the critica l density, actually incr eases smoothly from 
z = 4 —> 1.5 (D'Odorico et al. 2010) and maps well onto 
the, z < 1 values, measured with Hu bble Space T elescope 



(HST) ultraviolet spectra (Cooksey et al. |2010). Thus 



Ociv increases by, approximately, a factor ot tour over 
z rj 3 — > while dNcw/dX, the co- moving C IV line 
density, increases by a factor of two, roughly. 
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Early infrared spectroscopy, which probes the z > 5 
universe, first resulted in a co ntinuation of the roughly 
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constant ^cr y out to z k, 6 (Ryan- Weber et al. 2006 



Simcoe||2006 [). However, with limited sightlines, these 
studies were highly susceptible to cosmic variance, as 
shown by later IR surveys, which reported that £lciv 
actually dropped a t z ss 6 (Ryan- Weber et al. 2009 
Becker et al. 2009). The latest and largest IR survey 
showed that Vjciv drop s by, approximately , a factor of 
four over z w 4 — > 6 (Simcoe et al. 20111, and high- 
redshift quasars are continuing to be observed with the 
new Folded-port InfraRe d Echellete (FIRE) o n the Mag- 
ellan /Baade Telescope ( Simcoe et al. 2010). In addi- 
tion, Simcoe et al. ( 201ip found that all previous high- 
redshitt measurements overestimated flciv by ~ 30%, 
since their lower-resolution IR spectra led to the lower- 
redshift Mg II AA2796, 2803 being mis-identified as a 
strong, high-redshift C IV system. 
The recent low (UV) and high (IR) redshift publica- 



tions (Cooksey et al. 2010 
us to assess the state of t 
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ious studies have disparities in their definition of an ab- 
sorber, completeness corrections, sensitivity limits, and/ 
or adopted cosmology. Ideally, there should be a large, 
uniformly constructed, < z < 6 sample in order to 
evaluate the evolution of the C IV absorbers. We aim 
to produce an intermediate redshift catalog that is fairly 
comparable to the recent z < 1 and z > 5 catalogs. 

For our 1.5 < z < 4.5 survey, we use more than 
26,000 quasar spectra i n the Sloan Digital Sky Survey 
(SPSS; [York et al.|2000j) Data-Release 7 (PR 7) database 
( Abazajian et al.|A)09 " Schneider et al.|2010 ), thus mak- 
ing this study the largest C IV survey — both in path 
length and in number of absorbers — to date. Others 



et al. 2005 



have mined SPSS quasar spec tra for H I (e. g., Prochaska 



MM 
et al.l 



Pieri et al.||2010b| ), Mg II (e.g., Nestor et al 
ter et al.|20(J6D, Ca II AA3934, 3969 (e g 



Prodi 

2006[|Cherinka fe^chulte-Ladbeck|2011_ 
A, A1031, 1037 (e.g., [Frank et al.|2010[|Pieri et'al.|2010a[ ) 



and TTvT 



However, this is the first time that anyone has system- 
atically searched for C IV, likely because it is more diffi- 
cult and labor intensive, due to the increased amount of 
blending. This is the first in a series of papers on metals 
in SDSS quasar spectra, where we will assemble and ana- 
lyze self-consistent catalogs of Si IV AA1393, 1402, Mg II, 
and Ca II absorbers. These species are commonly stud- 
ied in quasar absorption-line spectroscopy and typically 
trace the gas closest to galaxies. 

We explain how we construct our C IV sample in [J2] 
and our completeness corrections in Sj3] The main re- 
sults are detailed in Q and the discussion and sum- 
mary are in §<|5jj6j We adopt the WMAP5 cosmology: 
Hp = 71.9kms" V Mp c, n M = 0.258, and 0\ = 0.742 
( |Komatsu et al.|2009 |. 



2. CONSTRUCTING THE C IV SAMPLE 
2.1. SDSS Data-Release 7 Quasars 
We began our survey wit h the 105,783 sightlin es in 



the SPSS DR7 QSO cat alog |Schneider et al.|poTol also 



see 



Richards et al. 2002 for the SUSS quasar selection) 
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Figure 1. Redshift and (S/N) distribution of sightlines. Using 
the left-hand axis, the 2D histogram shows the median signal-to- 
noise and zqso space of the analyzed 26,030 spectra. The black 
and gray histograms give the redshift distribution for all spectra 
and for the 10,861 with confirmed C IV absorbers, respectively 
(right-hand axis). 

9200 A and resolution varying from R = 1850 to 2200 

and the reduced spectra 



SPSS spectra have wavelength coverage 3820A < A < 



(or 162 km s 1 to 136 km s 1 
are binned to a log-linear scale of 69 km s - pixel - . We 
immediately excluded t he 6,214 broad-ab sorption line 
(BAL) QSOs detailed in |Shen et allpOlT ). 

We limited the sightlines to those with zqso > 1-7 and 
median signal-to-noise (S/N) > 4 pixel - in the wave- 
length range sensitive to C IV absorbers (see Figure 
mi. This range depended on the quasar redshift and the 
sPSS wavelength coverage. 

Previous, smaller surveys for C IV doublets have 
searched from the Lya A1215 to C IV emission wave- 
lengths. However, we removed a portion of this path 
length close to the Lya emission to avoid possible con- 
fusion between the C IV doublet and any O I A1302, 
Si II A1304 pair. These latter two have a wavelength 
separation similar enough to the C IV doublet that auto- 
mated search algorithms naturally bring in a large frac- 
tion of false positives. Since we had so much available 
path length available, we simply excluded a comfort- 
able region around the O I, Si II "forest" in addition 
to the Lya forest (i.e., A r > 1310 A in the restframe 
of the quasar). To exclude absorbers intrinsic to the 
QSO or affected by local ionization, clustering, and/ 
or enrichment, we set the upper wavelength bound to 
1548 A(H-2qso)(1 + (5wqso/c), 13 where we initially used 
Svqso — — 3000 km s -1 , but for the main analyses, the 
limit is Svqso = — 5000 km s -1 . 

After the redshift, (S/N), and coverage cuts, we had 
26,168 sightlines to search (see Table fl]). Later, we ex- 
cluded 138 sightlines as "visual BAL" CJSOs for a final 
count of 26,030 spectra to analyze. 

2.2. Automated Continuum Fitting 

We generated individual quasar continua for all ps 
26, 000 spectra in the sample using an algorithm which 

13 Velocity offsets are defined as 5v = c(z — z ro f)/(l + z ro f)- 
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Table 1 

Sightline Selection 



Number Description 

105783 SDSS DR7 QSO catalog {Schneider et al.|2010j 

99569 Excluding 6214 objects in BAL QSO catalog | |Shen et al.|2011[ ) 

48260 Covering max[1310A(l + ^ Q so), 3820 A] < A < min[1548 A(l + «qsq)(1 + Sv QSO /c), 9200A] a 

26168 With zqso > 1-7 and (S/N) > 4pixel _1 in above wavelength range b 

26030 Excluding 138 visual BAL QSOs 

10861 With confirmed C IV doublets 

a The 1310 A limit excludes the contaminating region in the OI, Sill "forest." The fiducial value for Svqso — 
-3000kms _1 . 

k The details of these specific sightlincs arc given in Table M 



Table 2 

Sightline Summary 



(1) 
QSO ID 



(2) 
R.A. 



(3) 
Decl. 



(4) (5) 
ZQSO (S/N> 
(pixel" 1 ) 



(6) 
/bal 



(7) 



(8) 
A/" can d 



(9) 
A/civ 



(10) 
<5^c,rv 



52235- 
52143- 
52203- 
52203- 
54389- 
51791- 
52143- 
52203- 
52902- 
51791- 
52203- 



0750-082 
0650-199 
0685-198 
0685-439 
2822-315 
0387-167 
0650-178 
0685-154 
1091-546 
0387-093 
0685-134 



00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:00: 
00:01: 



09.38 
09.42 
14.82 
15.47 
24.83 
39.00 
50.60 
53.17 
57.58 
58.22 
25.14 



+13:56:18.4 
-10:27:51.9 
-01:10:30.7 
+00:52:46.8 
+24:57:03.3 
-00:18:03.9 
-10:21:55.9 
-00:17:32.9 
+01:06:58.6 
-00:46:46.5 
+00:00:09.4 



2.2342 
1.8449 
1.8877 
1.8516 
3.2137 
2.1249 
2.6404 
2.7571 
2.5551 
1.8973 
1.9739 



5.94 
4.68 
5.16 
8.64 
9.21 
6.61 
8.68 
7.55 
11.83 
13.27 
5.32 



1.33 
0.83 
1.24 
2.16 
1.48 
1.75 
1.56 
1.15 
1.78 
0.84 
1.52 



028 



083 

084 
037 



Note. — Column 1 is the adopted QSO identifier from th e spectroscopic modifie d Julian date, plate, and fiber 
number. Columns 2 through 4 are from the DR7 QSO catalog l |Schneider et al.|2010| . Column 5 is the median S/N 
measured in the region searched for CIV. The binary BAL flag /bal in Column ti indicates which sightlincs were 
considered BALs by at least one author (4) and which were confirmed by the authors as BALs to exclude (8). Column 
7 is the maximum co-moving pathlcngth available in the sightline. Columns S and 9 give the number of candidate and 
confirmed CIV doublets, respectively. Column 10 is the pathlcngth blocked by the A/c iv doublets in the sightline. 
(This table is available in a machine-readable form in the online journal. A portion is shown here for guidance regarding 
its form and content.) 



combined principle-component analysis (PCA) fits, low- 
frequency b-spline correction, and automated outlier 
pixel and absorption-line exclusion. Since we are inter- 
ested in C IV lines only, we limited the fit to A r > 1230 A, 
redward of the Lya forest. We discuss details of the 
method below so interested users may replicate it; the 
codes are also publicly available in the XIDL software li- 
brary. 14 Figure [2] shows example fits for eight spectra of 
varying (S/N) level. 

First, every spectrum was iteratively fit using a 50- 
vector basis set of PCA "eigenspec tra" from the SDSS 
DR1 quasar sample ( Yip et al.|2004 ). The fit coefficients 
and their errors were dete rmined using the algorithms of 
Connolly & Szalay (1999 Equations 4-6), with spectral 
pixels weighted by their inverse variance when quanti- 
fying the goodness-of-fit. The formal per-pixel error of 
the PCA continuum fit was negligible compared to the 
spectral shot noise. 

A pixel exclusion mask was then generated upon each 
successive fit iteration. The initial mask consisted of the 
lowest 30% of pixels in successive 50-pixel windows, to 
filter strong absorption features and avoid biasing the 
continuum low on the critical initial fit. Then, for each 
subsequent iteration, a fresh mask was constructed which 
excluded pixels falling either less than — 2cr or more than 



Sec http : //ucolick. org/~xavier/IDL/ 



3<T from the working continuum. This mask was fur- 
ther modified to include the neighboring pixels of each 
deviant pixel below the working continuum and exclude 
regions with three consecutive pixels above. The latter 
correspond to systematically errant regions of the fit that 
should actually be included, so they were unmasked auto- 
matically. This cycle was repeated with successive PCA 
fits until: (a) the mask converged so successive fits did 
not change; (b) the unmasked flux and continuum fit had 
median difference less than 0.001%; or (c) ten iterations 
passed. 

For high-S/N spectra, low-frequency residuals in the 
PCA-normalized spectra become noticeable and affect 
our ability to search for weak C IV lines in the best data. 
We therefore performed a secondary correction to spectra 
with (S/N) > 8 pixel" 1 , re-fitting the PCA-normalized 
spectrum with an additional third-order b-spline, with 
breakpoint spacing of 25 pixels (or 1725 km s -1 ). As be- 
fore, we used inverse- variance weighting and clipped out- 
liers using the same thresholds. We refer below to this 
PCA plus b-spline fit as the "hybrid-continuum." 

The b-spline systematic error was estimated as the 
median difference between the b-spline and the PCA- 
normalized (masked) spectrum in bins matching the 
breakpoint spacing. This error was typically of order 
1% of the shot noise in the data; it exceeded the formal 
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Figure 2. Examples of automated continuum fits. Several QSO spectra (black) and 1-cr error (red) are shown with their automated 
continuum fits (green). On the left are the full spectra, and a zoom of the central wavelength range is given on the right. All spectra were 
fit for A r > 1230 A with SDSS eigenspectra. For (S/N) > 8 pixel -1 spectra, the eigen-fit was adjusted with a b-spline fit to the normalized 
spectrum. We surveyed the spectra for C IV absorbers between the (blue) vertical lines, which are outside the O I, Si II "forest" and more 
than 3000 km s" 1 blueward of the QSO. The gray brackets indicate real C IV doublets from the final sample. 



PCA fit error, and we added the two in quadrature to 
produce a total continuum error. 

The sigma clipping methods described above are ef- 
fective at identifying narrow features, instrumental arti- 
facts, or defects in the spectra. However, they are not 
always optimal for finding true absorption lines, which 
are always below the continuum; they can also be kine- 
matically complex and bias the fit. Since we have prior 
information about the characteristic wid ths o f cosmolog- 
ical absorption, we develop methods in §2.3| to machine- 
identify candidate absorption lines with these character- 
istics for the actual survey. We ran each PCA-normalized 
spectrum through a single pass of the absorption-line 
finder, using the resulting catalogs to generate a mas- 
ter absorption mask for final hybrid-continuum fitting. 

The full proce dure for automated feature finding is de- 
scribed in §2.3| briefly, we convolve both the data and 
error arrays with the instrumental response profile, and 
search for absorption features with signal-to-noise ratio 

> 3.5 per resolution element. For (S/N) < 16pixel _1 
spectra, we masked out all such features; in high-S/N 
cases, we also enforced that the unconvolved spectral pix- 
els at the line center must deviate from the continuum by 

> 5ct. Pixels falling within ±600 km s -1 of such features 
were added to the exclusion set. The mask consisting of 
automatically identified absorption lines was fed into the 
hybrid-continuum fit as a static mask. 

This procedure produced excellent continuum fits for 
the vast majority of spectra. When it failed, the most 
common reasons were: poorly measured zqsoj strong 
intrinsic absorption, broad-absorption lines, and fore- 
ground emission-line galaxy spectra superimposed on the 
quasar spectra. Instead of fixing these cases interactively, 
which would compromise our objective methodology and 
continuum error estimates, we chose to leave them in the 
sample and let the effects be accounted in our automated 
completeness and contamination tests (see Q. 

2.3. Automated Candidate Selection 



We limited our candidate search to C IV doublets 
where two automatically detected lines were within 
±150 km s^ 1 of the characteristic velocity separation 
<^civ = 498 km s -1 . The large uncertainty cut allowed 
heavily blended lines into our candidate list, but it also 
made the Mg II doublet our largest contaminant since 

^fMgii = 767 km s^ 1 . 

Our automated fea ture-finding algorithm is based on 
Prochter et al. (2006). The hybrid-normalized flux and 
error arrays were convolved with a Gaussian kernel with 
full-width at half-maximum equivalent to a resolution el- 
ement of the SDSS spectrograph, i.e., o~q — 1 pixel. The 
error array included the continuum fit error. The regions 
where the convolved signal-to-noise (S/N) conv is greater 
than or equal to 3.5 per resolution element were identi- 
fied as absorption features and saved to the sightline's 
line list. We masked out the few pixels (~ 5 A) around 
the strong skylines at 5579 A and 6302 A. 

Candidate C IV doublets were compiled by pairing 
automatically detected lines with the characteristic ve- 
locity separation <5tjciv ± 150 km s -1 and by identifying 
isolated, automatically detected lines that were broad 
enough to be a c andidate doublet by themselves. As 
discussed in |2.1| the search region is set so that the 
C IV candidate is redward of the O I, Si II "forest" 
(A r > 1310 A) and blueward of the QSO by 3000 kms" 1 . 
For any un-paired line, we re-ran the automated feature- 
finder with a 2.5er conv cut-off and searched for a partner 
candidate 1550 line. 

To identify isolated, broad lines that should be candi- 
date C IV doublets, we used the automated pr oced ure 
for measuring wavelength bounds (described in [2.4) to 
find lines with widths 5v\ iri > 1.58vcw- 

Ultimately, we identified 29,789 candidates, 6,346 of 
which were from the broad-line search. 

2.4. Measuring Absorber Properties 

The redshift, equivalent width, and column density of 
a line depend on the definition of its wavelength bounds. 
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The bounds were automatically defined by where the 
convolved signal-to-noise array stopped decreasing from 
the perspective of the automatically detected centroid, 
and the bounds were not allowed to exceed the mid- 
point between themselves and neighboring lines. Thus, 
for C IV lines, the inner wavelength bounds were not 
allowed to exceed the midpoint of the doublet. 

The redshift, equivalent width, and column density 
errors included the estimated error due to the hybrid- 
continuum. Rcdshifts were measured from the flux- 
weighted centroids. Equivalent widths were measured by 
simply summing the absorbed flux within the bounds. 

We used the apparent o ptical depth method (AODM; 
Savage & Sembach 1991 ) to estimate column densities. 
Systems that are saturated (discussed below) formally 
have column densities that are lower limits (binary flag 
/tv = 2). Column density measurements that are less 
than 3crjv are flagged as upper limits (/jv = 4). 

With doublets, we have two measurements of the col- 
umn density, and we use both to set New, as done in 
Cooksey et al. (2010). For doublets where both lines are 
measurements (not limits), New is the inverse- variance 
weighted mean of the two values. If one line is a mea- 
surement, New is set to its value. The column density 
limit is set to the inverse-variance weighted mean when 
both doublet limits are the same kind. When the doublet 
limits bracket a range, Nc w is set to the un-weighted av- 
erage with an error that reflects the range (/n = 8). In 
the remaining cases, New is set to the better measure- 
ment/limit and flagged /jy = 16. 

The AODM systematically overe stimates the tru e col- 
umn densities in low-S/N spectra ( Fox et aL]|2005 ). On 
the other hand, most doublets that can be detected in 
the low-resolution SDSS spectra are saturated, as evi- 
denced by the doublet ratios of near unity; therefore, 
the AOD column densities would formally be lower lim- 
its. We define the doublet to be unsaturated when the 
equivalent width ratio W r) i548/Wr,i550 > 2 — ctr, where 
ctr is the ratio error due to the uncertainties in the mea- 
sured equivalent widths. All other doublets are flagged 
as saturated with column densities that are lower limits; 
by this criterion, approximately 84% of the final catalog 
are saturated (see Table [3]). 

2.5. Interactive Doublet Verification 

The candidate C IV selection relied solely on the char- 
acteristic wavelength separation of the doublet. How- 
ever, true, unblended C IV doublets have a well-defined 
doublet ratio Wr,i548/Wr,i550 in the range of one to two 
in the saturated and unsaturated regimes, respectively. 
The resolved profiles of the doublet lines are very similar 
in the absence of blending. Most C IV absorbers are as- 
sociated with Lyq absorptio n, though "naked" systems 
may exist ( |Schaye et al. 2007| ) . There are frequently other 
metal lines associated with the C IV absorption, such as 
the Si IV and/or Mg II doublets. In addition, outside of 
the Lya (and O I, Si II) forest, line confusion and blend- 
ing are less severe, and the C IV doublet is one of the 
more common metal lines. 

We developed a graphical user interface to present all 
this information for each system to assist us in rating 
the 29,789 candidates. We leveraged our experience with 
C IV absorption systems and the varying amounts of 
information available, depending on redshift, to assign 
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Velocity (±900 km s 1 ) 

Figure 3. Example of C IV absorbers. One hundred doublet 
profiles are plotted from — OOOkms -1 to 900kms — 1 with respect 
to the center of the doublet and from —0.05 to 1.4 in normalized 
flux units. The flux is show in black and the 1-<t error in red. The 
horizontal dashed lines indicate the flux at unity (green) and zero 
(blue). 

each candidate a rating: 

0. definitely false — either can definitively identify the 
lines as other metal lines or as spurious bad pixels 
masquerading as absorption; 

likely false — though cannot name alternate metal- 
line identification, better data would probably con- 
firm this as not C IV absorption; 

likely true — though sparse supporting evidence 
(e.g., associated lines), better data would probably 
confirm this as C IV absorption; or 

definitely real — associated with other lines and/or 
shows clear correlations in line profiles. 
A single author was trusted to accurately assign ratings 
of or 3; these doublets were not viewed again. The 
more ambiguous cases (ratings 1 and 2) were reviewed 
by at least one additional author until consensus was 
achieved. We grouped the doublets into systems with 

o~Vzbs < 250kms~ 1 . 

Tho ugh we excluded the BAL QSOs from |Shen et ah] 
(20111, we found several sightlines with strong, self- 
blended, highly blueshifted C IV absorbers. Since we 
focused on the intergalactic absorption systems, we ex- 
cluded these 138 sightlines from further analysis and la- 
beled them "visual BAL" QSOs in Table [2] (binary flag 
/bal = 12). In addition, the extremely strong absorp- 
tion lines in these spectra corrupted the automated con- 
tinuum fitting algorithm. 

The final catalog includes all doublets with rating > 2. 
From the initial 29,789 candidates, 743 were in sight- 
lines excluded as visual BAL QSOs. In the remain- 
ing 29,046 candidates, we found 16,459 real C IV dou- 
blets (see Table [3]). A sample of doublets are shown 
in Figure [3] Ultimately, wc an alyzed the 14,772 with 
^QSO < -5000 km s" 1 (see 

We assessed the effects of blending and bad continua by 
inspecting 4500 doublets chosen at random. We flagged 
doublets where the 1548, 1550, or both lines were blended 
or otherwise problematic (e.g., the bounds were unrealis- 
tic) and where the continua ought to be adjusted locally 
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(e.g., due to strong absorption or neighboring emission 
lines). The majority of our analysis depended on W r ,i548j 
so it mattered most what fraction of the 1548 lines were 
troublesome. 

Over 90% of the random sample had "perfect" treat- 
ment of the 1548 line, and over 95% had acceptable con- 
tinuum fits. The problematic doublets showed some red- 
shift dependence due to the confused nature in the sky 
line region A G b s ^ 7500 A or 21543 > 3.8, where the frac- 
tion of problematic 1548 lines jumped from « 10% to 
« 20%. The stronger lines (Wr.i^s > 2 A) were more 
often blended. These systems tended to be self-blended 
and would generally be problematic to surveys using box- 
car summation for equivalent widths, such as ours. To 
maintain the largely automated, and hence objective and 
repeatable, nature of the survey, we decided to accept 
the imperfections and assess the strength of our results 
in light of them. 

3. COMPLETENESS TESTS 

We test our survey completeness by generating a li- 
brary of synthetic profiles, randomly distributing them in 
a subset of the sightlines, and tracking which simulated 
doublets were recovered. The goal was to populate a grid 
of doublet redshifts and rest equivalent widths with the 
fraction of recovered and accepted simulated doublets in 
each grid cell. We aimed to sample the full distribution 
of completeness limits so that we could leverage the large 
number of sightlines; for example, two spectra with 50% 
completeness in one (W r , 2154s) cell are equivalent to one 
sightli ne with 100% complete ness for the same parame- 



ters ( Matejek fc Simcoe][2012 ). In this section, we detail 
how we generated synthetic profiles, input and recover 
them, and measure the completeness, in addition to dis- 
cussing biases. 

3.1. Simulated Doublets 

We parameterized the simulated doublets to produce 
absorption lines that look like the diversity of profiles 
observed in the SDSS spectra. We created a Voigt profile 
library with a uniform distribution of equivalent widths, 
since we aimed to generate uniform errors on detection 
completeness as a function of W r ^4s- Each component 
had a column density that was sampled linearly from 
11 < logiVciv < 14 and a Doppler parameter randomly 
selected in the range 5kms _1 < b < 40kms _1 . For each 
system, one primary component was generated, and a 
random draw from a Poisson distribution, with mean of 
Mcompj determined the number of additional components 
for the system. The velocity offset for each component 
was drawn from a Gaussian distribution with standard 
deviation cr^f uss , clipped to keep all components within 
the maximum offset <5iw x relative to the system redshift. 



The bulk of our profiles were generated with /1 



or 7, cr 



Gauss _ 



OISS A 

comp 



100 kms~ 



and Sv„ 



= 300 km s" 



which were r oughly based on the com ponent properties 
measured by Boksenberg et al. (2003). However, to re- 



produce the diversity of profiles found in the SDSS data 
we had to modify the parameters (e.g., larger /U^^p and 
iu ma j.), and we limited the alternate parameters to the 
larger equivalent width r egime (e.g., W T > 1.2 A), where 
Boksenberg et al. ( 2003[ ) did not have many systems. 



We generated over 10 5 Voigt profiles. Then we ran- 
domly selected or duplicated (in extreme equivalent 
width regimes) profiles to uniformly sample 0.05 A < 
W t 1548 < 3.25 A in bins of 0.05 A, resulting in a library 
of 32,000 profiles. 

3.2. Monte Carlo Procedures 

We have two Monte Carlo completeness tests, referred 
to as basic and user. The basic test included all of the 
automated procedures, from continuum fitting to can- 
didate selection. It was applied to the largest fraction 
of sightlines (« 30%). We tested the user bias by visu- 
ally inspecting the candidates in sa 6% of the sightlines. 
The subsamples were set by requiring at least one sight- 
line in each bin of Azqso = 0.25 for zqso > 1-7 and 
A(S/N) = O^pixcr 1 for (S/N) > 4pixer 1 . For the 
basic test, with more sampled sightlines, we kept sam- 
pling from the sightlines in these Azqso and A (S/N) 
bins, so that the bins with more sightlines had more 
tested for completeness. Extra high-redshift sightlines 
were included in the user test in increase the statistics. 

We chose to input the simulated doublets into the 
actual spectra in order to sample realistic data. We 
"cleaned" each spectrum by removing a random 30% of 
the automatically detected absorption feat ures. The flux 
between the wavelength bounds (see {2 A) was replaced 
by continuum with a scatter drawn from neighboring pix- 
els, and the new error reflected those same pixels. By 
leaving the remaining absorption features in, we were 
able to measure the effects of blending and of absorption 
lines on the automated continuum fit and line-finding 
algorithm. 

We set a fiducial absorber redshift density dj\f civ jdz — 
5 to determine how many simulated doublets to input 
into the cleaned spectrum. This translated to typically 
injecting one to three doublets per loop, and the test it- 
erated until at least 1000 simulated doublets were input 
per sightline. We did not change the error array when in- 
jecting profiles, which slightly but systematically lowered 
the signal-to-noise ratio of the simulated doublets. 

For every sightline in the basic test, a random sam- 
ple of profiles were drawn from our library and assigned 
redshifts to fall within the limits of the current sightline. 
The spectrum was "cleaned" to start and every 50th loop 
thereafter. We injected one to t hree simulated doublets, 
fit the hybrid-continuum (see \2.1 1, and ran the auto- 
mated candidate identification algorithm (see £2.3). Any 



injected profile was flagged as recovered if the automat- 
ically identified candidate bounds spanned the observed 
wavelength of the simulated system. For any profile that 
was not recovered, we measured the flux-weighted red- 
shift, equivalent width, and AOD column density at the 
expected location. The input and re cove red information 
was stored for later processing (see E 3.3 ) . 

The user completeness test served a dual purpose. It 
tested the effects of human bias (e.g., we were less likely 
to accept real doublets at high redshift due to poor sky 
subtraction at > 7000 A) and our accepted false-positive 
rate (i.e., how often we rated non-C IV lines as C IV ab- 
sorption) . The steps were largely the same for one sight- 
line in the user test as in the basic, but we "cleaned" 
the profile every iteration and an author rated the au- 
tomatically detected candidates. However, we modified 
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Table 3 

C IV System Summary 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


QSO ID 


2QSO 


21548 


W r ,1548 


Wr,1550 


C(Wr,1548) 


log N c iv 


In 






(A) 


(A) 




(log(cm- 2 )) 





52143-0650-199 
52143-0650-178 



1.8449 
2.6404 



52902-1091-546 2.5551 



51791-0387-093 
52235-0750-550 
54389-2822-423 
54327-2630-423 
54452-2824-554 



52143-0650-561 
52991-1489-581 
52203-0685-567 



1.8973 
2.6383 
2.7668 
2.7251 
2.5566 



1.7593 
1.8536 
2.0921 



1.52755 
2.23461 
2.43947 
2.34348 
2.42902 
1.77015 
2.16695 
2.36223 
2.23946 
2.08349 
2.08788 
2.14618 
2.21361 
2.38185 
1.63718 
1.69764 
1.70301 
1.75102 



0.567 ± 
0.425 ± 
1.050 ± 
0.734 ± 
1.053 ± 
0.388 ± 
0.581 ± 
0.738 ± 
0.489 ± 
0.837 ± 
0.971 ± 
0.769 ± 
0.381 ± 
0.395 ± 
0.393 ± 
1.013 ± 
1.165 ± 
0.875 ± 



0.147 
0.059 
0.078 
0.202 
0.191 
0.097 
0.123 
0.087 
0.100 
0.040 
0.039 
0.041 
0.039 
0.039 
0.049 
0.090 
0.186 
0.121 



0.981 ± 
0.303 ± 
0.653 ± 
0.898 ± 
0.733 ± 
0.479 ± 
0.594 ± 
0.332 ± 
0.390 ± 
1.050 ± 
0.398 ± 
0.462 ± 
0.254 ± 
0.284 ± 
0.185 ± 
0.812 ± 
0.669 ± 
0.827 ± 



0.162 
0.058 
0.072 
0.198 
0.196 
0.098 
0.127 
0.088 
0.102 
0.042 
0.039 
0.037 
0.041 
0.039 
0.053 
0.088 
0.180 
0.122 



0.562 
0.363 
0.860 
0.713 
0.860 
0.305 
0.576 
0.716 
0.458 
0.781 
0.836 
0.738 
0.294 
0.317 
0.314 
0.849 
0.877 
0.799 



±0.206 
±0.092 
±0.023 
±0.192 
± 0.066 
±0.132 
±0.161 
± 0.073 
±0.151 
±0.025 
±0.014 
± 0.030 
±0.059 
±0.062 
± 0.075 
± 0.030 
± 0.038 
± 0.070 



> 14. 

> 14. 

> 14. 

> 14. 
?14. 

> 14. 

> 14. 
14. 

> 14. 

> 14. 
14. 

> 14. 

> 14. 

> 14. 
13. 

> 14. 
14. 

> 14. 



30 ±0. 
10 ±0. 

53 ±0. 
75 ±0. 
63 ±0. 
01 ±0. 
23±0. 
34 ±0. 
16±0. 
55 ±0. 
46 ±0. 
39 ±0. 
00 ±0. 
05 ±0. 
99 ±0. 
55 ±0. 
57 ±0. 

54 ±0. 



14 
06 
03 
14 
20 
13 
09 
06 
11 
02 
02 
03 
05 
05 
07 
04 
08 
07 



Note. — For each sightline (identified in Columns 1 and 2), every confirmed doublet is listed by the rcdshift 
of its CIV 1548 line (Column 3). The rest equivalent widths of the CIV lines are given in Columns 4 and 5. In 
Column 6, we give the completeness fraction for the doublet from the whole survey average. The C IV column 
density (Column 7) is the combined value from the AODM measurements in both lines. The binary column 
density flag Jn is composed of: 1 — good to analyze; 2 — lower limit; 4 — upper limit; 8 — un-weighted average; 
and 16 — default to line value with greater significance (Column 8). (This table is available in a machine-readable 
form in the online journal. A portion is shown here for guidance regarding its form and content.) 



the simulated doublets to have an unphysical character- 
istic wavelength separation of 4.8 A or 924 km s -1 so that 
any candidate was most likely either simulated or spuri- 
ous. The fake AA1547, 1552 profile library was truncated 
to W Fl i547 < 1 A. We did not simulate other absorption 
lines because, fundamentally, we are conducting a blind 
C IV survey. 

3.3. Completeness Correction and Unblocked 
Co-Moving Path Length 

We combined the completeness tests in bins of redshift 
and onto grids of equivalent width so that we correct 
any detected absorber based on its completeness frac- 
tion. The estimated completeness fraction for a given 
grid point (z s , W g ) is simply Accept /A/input, the fraction 
of input C IV doublets that are recovered automatically 
and accepted by the user, in the given grid cell. The 
two completeness tests separately measured the recover- 
ing and accepting effects so that the final completeness 
fraction is estimated by: 

C U w \ — ^cc(* g ,W 8 ) Kccc P t(z g ,W e ) 

l s ' gj ~ Amputfev^) A/; oc (z g ,wy ' {1) 

where the right-hand product can be thought of as 
CbasicCuscr- The completeness grid C(z g ,W g ) collapses 
to a curve C(W g ) in a fixed redshift bin. The basic com- 
pleteness uncertainty <Tc7 baaic is estimated by the Wilson 
score interval for a binomial distribution (Wilson 1927). 
This confidence interval estimator is well behaved for 
small A/input and/or for extreme completeness fractions. 

The statistics on the user completeness were naturally 
smaller than for the basic, so we fit C useY (W r ) with the 
following model: 



in two redshift bins. The dividing z = 2.97 corresponded 
to the beginning of the sky line region at « 7000 A, 
with its resulting decrease in C user due to confusion, and 
matched the start of the highest redshift bin. The fit un- 
certainties were estimated from Monte Carlo re-sampling 
of the equivalent-width errors, and the final error on 
C(W r ) was propagated from CT(7 baBic and <Jc ust!r - We ex- 
trapolated the fits to larger equivalent widths when we 
calculated the full completeness fraction. The user com- 



pleteness test is discussed in detail in i 3.4 below. 

We scaled C(z g , W g ) by the fraction of the total path 
length not obscured by doublets with greater equivalent 
widths, in the redshift bin. All C IV lines blocked < 2% 
of the total survey path length (see Table [2| . 

In actuality, A/Input and A/'rec contain the information 
only for the profiles and sightlines actually sampled in the 
basic completeness test. Any sightline that did not have 
a measurement for one or more quantities was accounted 
for with the average of the sightlines actually tested in 
the requested Azqso an d A(S/N) bin. 

The unblocked co-moving path length 15 for each grid 
cell is calculated by simply multiplying the completeness 
fraction by the total path length available in each redshift 
bin: 



AX(W g ) = AX(z g )C(z g ,W g ) 



(3) 



C useI (W T ) = C (l - e 



(2) 



In Figure[4j the black curves are the grid values AX(W g ) 
and errors. 

The unblocked co-moving path length to which our sur- 
vey is sensitive for any given detected equivalent width 

15 The co-moving path length is related to redshift as follows: 

X(z) = 2 V / n M (i±z) 3 ±rV(3Q M ). 
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Figure 4. Results from Monte-Carlo completeness tests. The completeness-corrected co-moving path length as a function of rest equivalent 
width AX(W Ty i54s) is shown for the 11 redshift bins used in the current study. The black curve and errors are the completeness curves 
(Equation pi , and the gray points and errors are the observations (Equation HI. The horizontal, dashed line traces the maximum path 
length available in the bin, and the vertical, dotted line indicates the equivalent width where we are 50% complete. 



Wi is interpolated from the grid of AX(W g ): 

AX(VFi)=interpol(AX(W/ g ), W g , Wi) (4) 
a AX =mterpo\(a AX{Wg) , W g , Wi) + 

(AX(Wi)- AX(Wi±a Wi )j . 

The error on AA(Wj) accounts for the error in our com- 
pleteness correction (a ax) and for the uncertainty in our 
equivalent width measurement (cwO- We plot each dou- 
blet's AX(W rj i548) in gray in Figure Bj Clearly, the un- 
certainty in equivalent width dominates the error in the 
completeness correction. 

In all redshift bins, the typical completeness fraction 
reached 50% by W Ty i^4s ~ 0.6 A, a value we frequently 
use as the minimum in subsequent analyses. 

3.4. Automated and Interactive Biases 

The completeness curves do not monotonically increase 
to 100% at large equivalent width (see Figure [4]) . In- 
stead, they roll off at W^^s ~ 2 A. This feature results 
from the incompleteness in the broad (self-blended) C IV 
search and, to a much lesser extent (< 5%), from broad 



profiles being over-fit by the continuum algorithm, which 
led them to being missed in the automated search. 

The user completeness test measured our ability to cor- 
rectly rate real, automatically detected doublets ("true 
positives" ) as well as the rate at which we include spuri- 
ous pairs of lines as true doublets ( "accepted false neg- 
atives"). We injected A/input = 5021 fake doublets with 
W T < 1 A and Sv QSO < -5000 km s" 1 , and A" rC c = 3070 
were automatically recovered. Of these, we correctly 
rated Accept = 2489. The left panels of Figure [5] show 
the trends of the user completeness with redshift, spec- 
trum (S/N), and equivalent width. 

Generally, C user decreases with increasing redshift and 
decreasing (S/N). Sky lines become numerous at A > 
7000 A, and poor sky subtraction leaves features in the 
spectrum that mimic absorption lines (see Figure pi. In 
the visual verification step, a larger fraction of real dou- 
blets are rejected due to the severe confusion. Lower 
signal-to-noise spectra induces the same effect. 

There were an additional 966 spurious candidates 
brought forth by the automated search, of which we in- 
correctly accepted 121. The accepted false-positive frac- 
tion peaks at ps 20% at z s» 3 and (S/N) pa 10. However, 




Figure 5. Biases of visual verification and trends of spurious detections. Left: We tested the "user bias" by rating fake AA1547, 1552 
doublets that were injected and automatically recovered as candidates. The four panels show the completeness of our ability to correctly 
rate true doublets, as a function of rcdshift, spectrum (S/N), and equivalent width in two redshift bins. The redshift cut (vertical line) 
was determined by the decrease in C uae r(^i547) due to the sky-line region at 7000 A and by the highest redshift SDSS bin. The (red) 
solid and dashed lines are the best-fit model (see Equation pi and its 1-u errors. Right: At the same time, we estimated the accepted 
false-positive fraction as functions of redshift, signal-to-noise, and equivalent width. These spurious pairs of lines were not injected but 
were automatically recovered. The lower panels show the best-fit model and errors to the accepted spurious fraction over W r ,i547- 

our acceptance fraction grows sharply at W r < 0.6 A and 
plateaus to 18% and 12% for z < 2.97 and > 2.97, re- 
spectively. The estimated co-moving line densities of ac- 
cepted false positives with W T > 0.6 A are dAf a f p /dX — 
0.022 and 0.020 for the low and high redshift bins, re- 
spectively. 

We detail how we applied corrections for the accepted 
false-positive rate in the next section. 

4. RESULTS 



The bulk of the analysis of the C IV sample was per- 
formed either on the whole dataset or in bins in rcdshift 
space. The bins were determined empirically from the 
$ v QSO < ~ 5000 km s -1 sample to have about 1500 dou- 
blets per bin, and the <5i>qso cut is explained below. 

4.1. Frequency Distribution 

The equivalent width frequency distribution /(Wr) is 
the number of detections A/obs(W r ) per rest equivalent 
width bin AW r per the total co-moving path length avail- 
able, in the given equivalent width bin, AX{W r ): 



maximum extent of the 1-er error ellipse on the likelihood 
surface (see Figure [7]) . All frequency distributions were 
fit over the range 0.6 A < Wm548 < max[W rj i548 + cr^/J. 
The results were not strongly dependent on the choice of 
the upper limit. The exponential model is a very good 
description of the data, and the best-fit parameters are 
given in Table |4j 

The best-fit parameters show smooth redshift evolu- 
tion with respect to the constant dAfciv /dX curves in 
Figure [7] The co-moving line density is simply the in- 
tegral of the frequency distribution from some limiting 
equivalent width Wu m to infinity, and substituting in our 
exponential model, we see: 



dtft 



civ 



dX 



-k 



aWii, 



(7) 



lit 



f(W t ) = 



AW r AX{W r ) 



(5) 



/(Wr) - 



■■f(Wr) 



\AX(W r 



The error on Af b s is estimated from a Poisson distribu- 
tion if A/" bs < 120 and from a Gaussian approximation 
if A/"obs > 120. For AX(W T ) and a ax, we used Equation 
[4] with the center of the equivalent width bin being W 
and o- Wl = 0.5AW r . 

We used t he maximum likelihood analysis of |Cooks"ey| 
et al. (2010) to fit /(W r ) with an exponential: 



By fixing dA/civ I 'dX , we can solve for the required nor- 
malization k for any given a. The 1-a ellipses are elon- 
gated in the direction of the constant dAfciv /dX curves, 
and redshift evolution can be seen by tracking system- 
atic change perpendicular to these curves. The lowest 
five redshift bins (1.46 < z < 1.96) fall along roughly the 
same constant dAfciv /dX curve. The next three highest 
redshift bins (1.96 < z < 2.51) have slightly smaller line 
densities. Then there is almost a factor of two drop over 
the highest two redshift bins (2.51 < z < 4.55). 

Since the accepted fals e-po sitive rate is essentially con- 
stant at W r > 0.6 A (see { 3.4 ), the frequency distribution 



of accepted false positives is a scaled-down version of the 
measured frequency distribution. Therefore, we scale the 
original /o(Wr) as follows: 



/(W r ) = 1 - 



f{W v ) = ke aW < (6) 

(see Figure p| . The normalization k and scale a were 
simultaneously fit, and the errors were estimated by the 



dAf afp /dX\ 

Jdx) f ° m 



dAf ( 



civ I 



(8) 



and propagate the errors. This results in a decrease of 
« 6% to 12%, depending on the redshift. For the expo- 
nential fits, we scale the best-fit normalization fco in a 




Figure 6. Equivalent width frequency distributions. The maximum likelihood fits of an exponential function are the dashed (red) lines, 
for each redshift bin, and the solid (gray) line, for the full sample. There is little evolution with redshift. The observations have been 
completeness corrected, and the redshift-specific 50% completeness limits are the vertical dotted lines. 



similar fashion: 



dM aip /dX 



in 



(9) 



but the denominator is the integrated line density from 
the exponential model (Equation^]). We report the prop- 
agated errors in Tableland in the text. However, since 
we cannot compute the change in the likelihood surface, 
we only shift the ellipses in Figures [7] and [8] 

4.1.1. Effect of Blending and "Intrinsic" Absorbers 

As mentioned in §2.5| blended profiles are an issue for 
this survey, which relies on automated procedures and 
boxcar summation to measure equivalent widths. We 
ran 1000 Monte Carlo simulations to estimate the net 
effect of blending. There were over 10 6 (17%) simu- 
lated doublets recovered in the completeness tests that 
had measured equivalent widths more than 3cr larger 
than the input value, indicative of blending. We mea- 
sured the median ((Wbiend) = 0.2 A) and standard de- 
viation (criy bkmd = 0.5 A) of the distribution, where 
Wblend = Wr.rec - (Wr,in P ut + ^W t , re J- For each real- 
ization, we tested the worst-case scenario by decreasing 
the equivalent width of a random 20% of the absorbers, 



which was the largest blended fraction estimated in j 2.5 
The magnitude of the decrease was drawn randomly from 
the half of a Gaussian distribution below its mean, set 
to (Wbiend), with standard deviation crw bl<md . Then, the 
new sample was fit with an exponential f(W r ) model. 
The median best-fit parameters from these Monte Carlo 
simulations were in very good agreement with those in 
Table |4j and the standard deviation of the Monte Carlo 
results were ten-times smaller than the quoted uncertain- 
ties. These simulations assessed the worst-case scenario, 
by ta king the largest estimated fraction of blended lines 
(see [2.5| and only decreasing the equivalent widths (as 
opposed to also re-sampling the the other 80% of the 
doublets). 

In addition, the completeness corrections account for 
blending since we compiled the curves in Figure[4]with re- 
spect to the measured — as opposed to input — equivalent 
widths. There were other absorption lines in the spectra, 
and the randomly placed simulated profiles were blended 
at a realistic rate, as seen by the ag reem ent of the blended 
fraction from our visual estimate ({ 2.5 ) and from the 10 6 
simulated doublets described above. Thus, our results 
are robust to the effects of blending. 

We tested the effect of the Svqso cut on our re- 
sults by increasing <5i>qso and re-fitting the full red- 
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Table 4 

C IV Results Summary 



(1) 


(2) 


(3) 


(4) 


(5) (6) (7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


<*> 


■^lim 


A/"obs 




W 50% dMciv/dz dMcw/dX 

(A) 


S^C IV 


Afet 


k 


a 












(xlO- 8 ) 




(A- 1 ) 


(A- 1 ) 



1.96274 
1.55687 
1.65971 
1.73999 
1.82360 
1.91460 
2.01778 
2.15320 
2.35608 
2.72298 
3.25860 



[1.46623, 
[1.46623, 
[1.61004, 
[1.70035, 
[1.78007, 
[1.87000, 
[1.96002, 
[2.08009, 
[2.24015, 
[2.51028, 
[2.97005, 



4.54334] 
1.60986] 
1.69993] 
1.78000] 
1.86985] 
1.95998] 
2.07997] 
2.23997] 
2.50914] 
2.96976] 
4.54334] 



14772 
1428 
1473 
1449 
1568 
1424 
1523 
1490 
1450 
1503 
1464 



38624 
3962 
3420 
2987 
3067 
2665 
2983 
2878 
3341 
5673 
7632 



0.52 
0.63 
0.57 
0.53 
0.51 
0.49 
0.46 
0.45 
0.50 
0.61 
0.59 



0.92 
0.95 
1.00 
1.08 
1.09 
1.12 
1.04 
1.04 



+0.02 
0.01 
0.03 
0.03 
+0.04 
0.03 
+0.04 
0.04 
+0.04 
0.04 
+0.04 
0.04 
+0.04 
0.04 
+0.04 



-0.04 

96+ 04 
u - yD -0.03 
n OO+0.03 
U - Scl -0.03 
59+ 02 



n 97c; +o.oo4 

u ' z '°-0.004 
n oon+0.011 
u.ooz_ Q Q10 

336+ 012 
u.ooo_ 011 

U.OOD_ 012 
n ocri+0.012 
U " )oi -0.011 
n n^+0.013 
u.ooo_ 012 

n o 99 +0.012 
U ' J "-0.011 
n o 19 +0.012 
u - olz -0 .011 

276+ - 011 
u - z ' D -o.oio 

n 991 +0.009 
"■^ -0.008 
14 r+0.006 
u - 1 ^°-0.005 



1.71 
2.18 
2.10 
.17 



.15^ 



2 
2 
2 
1 

1.92 
1.61 
1.44 
0.87 



+0.20 
0.20 
0.3!) 
0.39 
+0.32 
0.32 
0.28 

0.28 
1-1+0.23 
11 0.24 
0.23 
0.23 
+0.17 
-0.17 
+0.17 
0.17 
+0.18 
0.18 
+0.22 
0.22 
+0.13 
0.13 



8918 
1025 
956 
902 
932 
824 
852 
807 
775 
969 
876 



3 49+0.31 

o 87 +0.94 
-0.82 

77 +0.95 
°" ' ' -0.82 

4 46+ 1 - 24 
^■ 4D -1.06 

1 40+1.20 
^■* 3 -1.04 

B.l8ii;S 

,+1.21 
-1.04 

4 2fi +1 ' 62 
*' ZD -1.37 
i+0.74 
-0.65 
,+0.75 
0.65 



4.0V 



2.48] 
1.82"' 



-2 62+ - 04 
z - uz -0.04 

9 cn+O.U 

z ' oz -0.12 

-2 48+ - 12 
z -*°-0.13 
1+0.14 

-0.14 

-2 63+ 13 
z - uo -0.13 
,+0.14 
-0.14 
9 S7 +0.15 
- Z -S<_0.15 
,+0.14 
-0.14 
_ 9 S9 +0.15 
z - oz -0.16 
-,+0.12 
-0.12 
,+0.13 
0.13 



-2.6r 



-2.63"* 



-2.63"* 



-2.49Z 
-2.61 H 



0.577 
1.560 
0.876 
1.672 
1.291 
1.832 
2.838 
3.267 
1.296 
0.837 
1.548 



Note. — Summary of the most common redshift bins and data used for the various analyses. Columns 1-2 give the median, minimum, 
and maximum rcdshifts for the observed number of doublets (Column 3), and the maximum co-moving pathlcngth in the redshift bin is given 
in Column 4. The 50% completeness limit from the Monte Carlo tests is in Columns 5. The redshift and co-moving absorber line densities for 
W T > 0.6 A are in Columns 6—7. In Column 8, the Ociv from summing the mass in the W T > 0.6 A absorbers is a lower limit, since the majority of 
absorbers arc saturated. The frequency distribution was fit with an exponential /(VK r ) — feexp(ctWr) for A/fit absorbers with W T > 0.6 A (Column 
9), and the best-fit parameters are given in Columns 10—11. The reduced x 2 from the best fit and f(W T ) (in bins with ~ 100 doublets each) is 
given in Column 12. 
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Figure 7. Best-fit f(W r ) parameters and errors. We fit the fre- 
quency distribution with an exponential function (see Equation^. 
The best-fit normalization k and scale a and the l-cr error ellipses 
are plotted for the ten small redshift bins (numbered points); the 1, 
2, and 3-c contours are shown for the fit to the full sample (black 
ellipses). The best-fit parameters smoothly evolve with redshift, 
as seen by comparing the ellipses with the constant dMcw/dX 
curves (gray, dashed lines). 

shift sample (see Figure [8]) . Systems close to the back- 
ground QSO could potentially be high velocity, intrin- 
sic absorbers and/or affected by quasar clustering, en- 
richment, and/or local ionization. The general trend for 
Svqso < —4000 km s~ is decreasing normalization (i.e., 
k) with nearly constant shape (a). The decrease in k is 
due to the number of absorbers decreasing faster than the 
path length up to at least Svqso = — 15000 km s -1 . For 
example, with Svqso — — 10000 km s~\ we have 17.5% 
less path length but 25% fewer W T > 0.6 A doublets. 
This trend is partially ex plained by quasars res iding in 
dense environments (e.g. 



Prochaska et al. 



rp011| ), where 
the chance of intersecting a metal-enriched galaxy halo 



Figure 8. Best-fit f(W T ) parameters for different Svqqo values. 
We tested the effects of the (Svqso cu t on the fit to the full C IV 
sample. The contours are 1 and 3(7, and the gray, dashed lines are 
constant dAfciv/dX curves. We adopt (Sdqso = — 5000 km s -1 for 
the main analyses. 

is increased. 

However, this effect would not dominate out to 
Svqso = — lOOOOkms -1 . The lack of "convergence" 
at very large Svqso may be partially due to po orly 
measured quasar redshifts. |Hewett fc W ild (20 10| re- 
measured the m ajority of the redshitts for the |Schnei-| 
der et al. (2010) DR 7 QSOs, with better au tomated rou- 
tines. Adopting the |Hewett fc Wild| ( |2010| ) redshifts af- 
fects over 80% of the total 16,459 C IV systems. The 
majority of the quasar redshifts increased, with a median 



(Sv c 



381 km s 1 , standard deviation 764 km s 1 , and 



maximal extent — 4200kms _1 < 5v em < +5700 km s -1 . 
In addition, ther e are redshift-dependen t fluctuations in 
Sv em . Using the Hewett & Wild (20101 redshifts would 
affect our analysis; however, the change in the total sam- 
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Cooksey et al. 



pie size, given a Svqso cut, is less than 2%. Given this 
small fraction and the incompleteness of the new red- 
shifts, we chose to continue with the Schneider et al. 
(2010) redshifts, but we adopt Svqso = —5000 km s~ 
tor the bulk of our analyses which reduces the sample to 
14,772 doublets. 
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Figure 9. Co-moving C IV line density evolution. The number 
of absorbers per co-moving path length increases steadily and con- 
sistently from 2i548 = 4.5 — > ~ 1.74. As expected from the nearly 
unchanging nature in the shape of f(W r ) (see Figure pi, there is 
little dependence on W r ii m , noting that we are typically 50% at 
W T « 0.6 A. 



4.2. C IV Absorber Line Density 

We directly measured the absorber line density for dou- 
blets with W T > Wy lm as follows: 



dN< 



CIV 



dX 



(Wr > Wa 



limj 



T dM/dX 



AX(z) 
AX(z) 



(10) 



The completeness-corrected number of absorbers in any 
given bin is the completeness-weighted sum of the ob- 
served absorbers in the bin: 



WiKWr+O.SAWr 



AA c (w r ; 



E 



i 



Wi<W r +0.5AW r 



°M C = a lf ohs + 



C(Wi) 

r +0.5Aff r / x . 

n kaw V V 1 ' / 



(11) 



Wi>lV r -0.5AIV, 

Again, N \>& is the contribution of the actual observed 
number of absorbers in the given W r bin, and the 
completeness-corrected number A/"c > A/" bs- The error 
er/Voba is estimated from a Poisson distribution if Af bs < 
120 and from a Gaussian approximation if Af bs > 120. 

We subtract dN a f p /dX from all quoted A A/~ciy /dX val- 
ues, in the appropriate redshift bins (see [3.4) and add 
the errors in quadrature. For W T > 0.6 A, dM a i v /dX — 



-0.003 



0.020; 



0.020^ (z > 2.97). 



(1.4 < z < 4.6); 0.022+g;^! (z < 2.97); and 



We present dAfcw/dX for different equivalent width 
limits in Figure [9] The line density shows little dif- 
ferential evolution based on Wrj; m , as expected from 
the consistent shape of the frequency distributions over 
time (Figure [6]). For W r > 0.6 A, the inverse- variance 
weighted average <W C xv/dX = 0.350 ± 0.005 for 1.46 < 
z < 1.96 (or the lowest five redshift bins). This av- 
erage is 2.37 ± 0.09 times larger than the highest red- 
shift d/Vciv /dX at z = 3.76. Though the magnitude 
of the increase is modest, the detection is a > 20er 
result. Thus, the line density grows consistently and 
smoothly from z 1548 = 4.5 — > ~ 1.74, then plateaus at 
dNcw/dX « 0.34 until z = 1.46, as expected from the 
best-fit f(W T ) parameters (see Figure J7|. 

There is a known bias in the SDSa quasar color se- 
lection that leads to an excess of Lyman-limit systems 



at 3 < -zqso 5, 3.5 (Prochaska et al. 2010 Worseck & 



Prochaska 2011), which likely increases the incidence ot 



strong mefai-hne absorption systems at these redshifts. 
Therefore, there is potential for the decrease in the high- 
est SDSS bin to be even larger, but the effect of the color 
selection on the C IV sample is beyond the scope of this 
paper. 

4.3. Comparison with Previous Results 

To fairly compare to other surveys, we applied various 
W r cuts to our complete sample, chosen to match the 
corresponding cuts of prior surveys. 

To compare f(W r ), we converted p ublished column 
dens ity frequency distr i bution s jits from Songaila (2001) 



and D'Odorico et al 
quency distributions 



(2010) to equivalent widt 



fre 



ay assuming the linear curve of 
growth and mapping /(Aciv) directly to f(W T 154s): 
f(W T ) = f(N clw )dN cl y/dW r . We adjusted for cos- 
mology as necessary. A single-component cloud becomes 
saturated and nonlinear at logAciv ~ 14, which trans- 
lates to Wr,i548 ~ 0.6 A, where we typically are « 50% 
complete. 

While we have good statistics on the rare, strong ab- 
sorbers, we suff er f rom incompleteness at W r < 0.6 A. 



Songaila ( 2001 ) and D'Odorico et al. ( 2010 ) were smaller, 
higher-resolution, higher-S/JN studies, and they were 
complete to very low W T but suffered from sample vari- 
ance at larger equivalent widths. However, Figure 10 
shows that our results are consistent with each other 
in the overlap region. O ur best-fit parameters to the 
exponent i al f(W r ) in the D'Odorico et al. (2010) and 
Songaila| (2001 ) redshift bins were, respectively: a = 



-2 fi^+ - 04 A" 

^■OO_ Q 4 A 



1 and k = 3.72±g;|^ A 1 for 1.6 < z < 3.6 



with 7884 systems; and a 
2 29+ - 80 A" 1 

z - za -0.69 A 



-2.58; 



-0.13 X" 

-0.13 A 



and k 



for 2.9 < z < 3.54 with 878 systems 
We compared the SDSS red shift and co-moving line 



densities with the literature ( Steidel 



Tytler|[T998| [Cooksey et al. 



1990 
dorico 



Barlow & 



'Odorico et al. r 



Simcoe et al.|2()l ip. For |D'6dorico et al.| ( |2010| 7we usee 
their best-ht /(iVc iv) values to calculate d/Vc iv / dX and 
dM cw/dz = (dAfmy/dX) (dX/ dz) for l ogJVc Tv > 14. 



For |Barlow fc Tytler| ( |1998[ ) and |Steidel| ( |1990[ ) , we esti- 
mated dM cw/dX = {dTJcw/dz^dX/dzT^ The [SteT- 



del ( 1990 ) results match ours well for the range where the 
author was fairly complete but with w 2 0% uncertainty, 
compared to our w 2% errors (see Figure 11 ). Extending 
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Figure 10. Comparing f(W T ) to smaller, higher resolution studies. We show the best power-law fit from left: 1.6 < z < 3.6 (ct p i = 
-1.71 ± 0.07; p'Odorico et al.||2010} and right: 2.9 < z < 3.54 (ct pl = -1.8 ± 0.1; |Songaila|[200T| l as the (blue) dashed lines. Our best 
exponential fit for each redshift bin is shown on top of the observed f(W T ) (in bins of « 500 and 100, respectively). The vertical 50%- 
complete line also happens to be around where a single C IV component saturates. The gray plus-signs indicate the equivalent widths 
where the labeled column densities fall. 
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Figure 11. Reds hift (left) and co-moving (right) C IV absor ber line densitie s . We compare appl es-to-apples dNr.iy/dz with W T > 0.4 A 
with five surv eys: |Barlow &: Tytler] j!998| light gray squa re), [Cooksey et al,| ( |2bl0| green cross), [D'Odorico et al.j l |2010| black triangle), 
|Steidel| ( |1990[ gray diamonds), and |simcoe et al.| ( |2011| red crosses). Tne co-moving line density steadily increases trom high to low 
redshift. The growth in dAfciv/dX indicates a consistent increase with time in the product of the co-moving number density and the 
physical cross-section of the absorbing clouds. 

the redshift coverage by including th e lo w-redshift mea- 
surem ents of |Barlow fc Tytler ( 1998| ) and |Cooksey et ah 
and the high-redshift values from |Simcoe et al. 



2010 



2011 



we see that dNciv/dX has steadily increased 



volume density of absorbing clouds n com and their phys- 
ical cross-section <T p h y s : 



by, roughly, a factor of ten from zi^g = 6 — > 0. 

5. DISCUSSION 
5.1. C IV Evolution 
We show the best measurements of dAfcw /dX 



dAfcw 
dX 



H 



7 phys- (12) 



increase m 



for 



Wr > 0.6 A and < z < 6.5 in the top panel of Figure 



12 The co-moving line density relates to the co-moving 



16 For z < 1, we summed the e ight W T > 0.6 A systems in the 
HST sample |Cooksey et al.|2010 [l to get (W C iv/dX = 0.63l°". 
For z > 4, we have dM C lv/d,X = 0.212 ± 0.092 for z = 4.94 (2 



Thus, the roughly order-of-magnitude 
dAfcw /dX from high-to-low redshift means the product 
^comCphys has increased by a factor of « 10. We know 
the metallicity of the universe has steadily increased over 
cosmic time, so the possible number of C IV-absorbing 
clouds (i.e., n com ) has likely increased. At least some 
C IV absorption traces g alaxy halos at low (i Ch en et al. 
2001[ ) and high redshift ( |Adelberger et al.||2005| |Martin 



systems) and 0.026 ± 0.027 for z ■■ 
l |Simcoe et al.|2011^ . 



5.66 (1) from the FIRE sample 
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Cooksey et al. 



et al.|2010[ [Steidel et al.|2010[ ). Since galaxies (and their 
halos) have likely grown over cosmic time, increases in 



both (Tphys and n, 
i n dNc Yv/dX 



appear to contribute to the increase 



Adelberger et ah] p005| ) 
logiVciv > 14 absorbers arise within 



reported that almost all 
80 kpc of Lyman- 
break galaxies (LBGs) at 2 < z < 3. The evidence in- 
cluded: individual strong absorber-LBG pairs; strong ab- 
sorption in stacked spectra of close background galaxies, 
shifted to the rest-frame of the foreground galaxies; and 
similar LBG-C IV cross-correlation and LBG autocorre- 
lation functions, suggesting they have the same spatial 
distribution. 



Steidel et al. (2010) increased the LBG-LBG pairs for 
stacking analysis, and they measured average W r ,154s = 
0.13 ± 0.05 A and 1.18 ± 0.15 A at distances b = 63 kpc 
and 103 kpc, respectively. Interpolating between these 
two measurements, an average W r 154s = 0.6 A s ystem 
would reside at b « 85 kpc 
LBG sample went as faint as 



The iSteidel et 
T3T7 



al. (2010) 



ass uming the 

luminosity function of Reddy & Steidel ( 2009 ) , and since 
they used stacks of galaxy spectra, their analysis included 
the effects of partial covering fractions. 

Since the W T > 0.6 A C IV absorbers in our sample 
have log-/Vciv ^ 14, we used the co-moving number 
density of UV-selected galaxies, to estimate the typi- 
cal galaxy-C IV cross-section over time with Equation 
>tcom.uy with the UV luminosity func 



12 We measured 



"flons from Oesch et al. 



andlBouwens et al. 1200 



~aTT(| 2~ 



2010) 



Reddy fc Steidelj ( [2009 1 
1 covered several smaller 



whic 

redshift bins spanning 0.5 < z < 2, 1.9 < z < 3.4, and 
3-8 < z < 5.9, respectively. For each luminosity func- 
tion, we integrated the best-fit Schechter function down 
to 0.5 L* or sa 0.75 mag fainter than the published M^ v , 
which ranged between ss —19 and —21 mag. We esti- 
mated the n comj uv errors with Monte Carlo simulations. 

The UV-selected galaxy number density, n conii uv, may 
increase by a factor of two to three from z w 6 — >■ 0, but 



the uncertainties are large (see Figure 12 middle panel) 



12 



we 



Applying n conii uv and dNciw/dX to - Equation 
estimated (Tphys, the physical galaxy-C IV cross-section 
(Figure|l2j lower panel). For some redshift bins, multiple 
luminosity functions could be used in conjunction with 
our d/Vcrv /dX , and we show all resulting values but 
highlight the preferred values with filled symbols. 

The galaxy-C IV cross-section shows no evolution over 
the SDSS redshift range within the errors, which are 
dominated by the 20%-60% uncertainties in n com jjv- 
Assuming the cross-section is due to a spherical halo 
that projects with 100% C IV covering fraction, the 
halo radius would be i? p h ys = \/<7phys/7r ~ 50 kpc for 



1.5 



< 4.5. 



This distance agrees with Adelberger 
et al. ( 2005 1 and allows for e.g., a non- unity covering frac- 
tion or limiting to brighter UV-selected galaxies, both of 
which would increase i? p hy S - 

There is an approximately 10-fold increase in <T p hy S 
from z w 6 — > when we inc l ude t he d Ncw/dX mea- 



surem ents from Simcoe et al. (20111 and Cooksey et al. 
(2010), respectively. This increase in the galaxy-C IV 



cross-section is comparable to that of dNciv/dX over 
the same redshift range. Thus, if the W T > 0.6 A C IV 
absorbers were only tracing galaxy halos, the redshift 
evolution of dAfciv/dX could be solely due to the halos 



filling up with triply-ionized carbon, through some com- 
bination of physical growth, increased metallicity, and/ 
or evolution in the ionizing background. 

However, the uncertainty in the z < 1 cross-section 
estimate is large and consistent with an increase of only 
two to three compared with the z ss 6 value. In this case, 
the roughly order-of-magnitude increase in dMcw/dX 
could equally be due to modest increases in n comi uv and 

Cphys ■ 

We emphasize that the preceding discussion assumes 
all W T > 0.6 C IV systems at < z < 6 are only found in 
the circum-galactic media of UV-selected galaxies. Any 
intergalactic C IV contribution to dNcw/dX at any red- 
shift would cause us to overestimate the galaxy-C IV 
cross-section. 

Evolution in the ultraviolet background (UVB ) possi- 
bly ex plains the seeming spike in cr p h ys at z « 5. [Simcoe 



(2011 1 showed that C IV is a preferred, if not dominant, 
transition of carbon at z = 4.3 compared to lower r ed- 



shift, assuming the UVB of Haardt fc Madau 



Faucher-Giguere et al. ( 2009 ). Though the redshift win 



(2001) or 



dow of this effect is small (Az ~ 0.5), the increase in 
C IV-ionizing photons would increase the C IV cross- 
section of UV-selected galaxies and may explain the spike 
in d/Vcrv /dX and <x p h ys - In general, triply-ionized car- 
bon increasingly becomes a disfavored transition of car- 
bon with de creasing redshift, given a s t andard mode l 
of the UVB ( |Oppenheimer fc Dave|[2008| |Simcoe|[20lT ). 
The UVB likely dominates the ionization flux at the ex- 
pected tens of kiloparsecs from the host galaxies. Thus, 
the increasing number of C IV absorbers towards lower 
redshift indicates an increasing enrichment of the gas, 
to more-than-compensate for the decreasingly favored 
triply-ionized state. 

C IV a bsorption could als o be a galactic wind signa- 
ture (see Steidel et al. 2010). However, the dAfcw/dX 
redshift evolution does not mimic the evolution of 
dAfyigu/dX for W ri 2796 > lA absorbers, and these 
strong Mg II doublets are often analyzed as wind tracers 
(iWeiner et al.| |2009[ [Rubin et al. | |2010[ i Bordoloi et al 



201H but also see |Gauthier et al.||2010| and |Kacprzak 



et al.|20lT ). Strong Mg II absorbers evolve strong ly with" 



redshift ~5y increasing with incr easing redshift (Nestor 
et al.||2005| |Prochter et"al]|2006| , peaking 
then decreasing at higher redshifts (?). Tl 



at z « 3, and 
(?). The latter study 
showed that the evolutio n in dNyieiil dX trac ks the cos- 
mic star-formation rate (jBouwens et al.||2010 ), using the 



scaling relation of |lVlenard et al. l 201ip , reinforcing the 
idea that strong Mg 11 absorbers arise m galactic winds. 

We see no decrease at lower redshifts for strong C IV 
absorbers (see Figure |9| . If strong C IV and Mg II ab- 
sorbers were both tracing winds, then, to account for 
their different redshift evolution, either: they must probe 
outflows in different ways; or strong C IV absorption 
traces an additional medium, such as halo gas, that con- 
tributes significantly to its d/Vciv / 1 dX evolution. 

5.2. C IV Mass Density 
The C IV doublet redshifts into optical passbands for 



1.5 < 



< 



5.5 and has been o bserved extensively with 



large, ground-based telescopes (Songaila 2001 Boksen 



berg et aLl|20031 ISchaye et al.||2003| |Scannapieco et"aT 
al."" 



2006| |D'Odorico et al.||2010l)/ ~^urveys of ultraviolet 
quasar spectra taken witn HST cover the z < 1 C IV sys- 
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Figure 12. Evolving nature of C IV absorbers. Top: By com- 



bining the current SPSS v alues with those from |Cooksey et al. | 
| |2010| l and |Simcoe et al.| l |2011| , we show that dLNciv/d-X lor 
Wr,i548 > 0.6 A has increased 10-fold, roughly, from z = 6 — > 0. 
Middle: We calculated the co-moving number density of UV- 
selected galaxies by in te grating the U V luminosity functions from 
Bouwens et al.| \2007\ , |Oesch et aT] \2010\ , and |Reddy fc Stei-| 
del| j20Ulj| , down to 0*5 L~ . UV-selected galaxies are bright, star- 
lorming galaxies, such as the LBGs that [Adelbergcr et alij | |2005[ l 
showed likely host most log Nq jy > 14 absorbers at 2 ^ z ^ 6. 
Bottom: Assuming that all 14^,1548 > 0.6 A doublets arc in UV- 
selected galaxy halos, we estimated the gala xy- C IV cross-section 
as a function of redshift with Equation |12| The inferred cross- 
sections from the SPSS 1.5 < z < 4.5 sampl e are consistent with 
the area that |Adelberger et al.| (2005| and |Steidel et~aT"| | |2010[ l 
measured. 



recently, improved infr ared spec trographs have pushed 



terns (|Danforth fc Shull|2008||Cooksey et al.|2010[). More 

i ha 



et al. 


2009 e 


have 


focused 



Simcoe et a 



6 (Ryan- Weber et al. 2009 Becker 



Typically these studies 
ficiv, the C IV mass 
density relative to the closure density. The first moment 
of the column density frequency distribution /(-/Vcny) is 
related to ficiv as follows: 



H m c 



civ 



c Pc,0 Jn, 



AL 



f(N C iv)N cw dN cw , (13) 



where Hq = 71.9 km s~ Mpc~ is the Hubble constant 
today; mc = 2 x 10 -23 g is the mass of the carbon 
atom; c is the speed of light; and the critical density 



pc,o = iHf^G)- 1 = 9.77 x lO- 30 gem" 3 for our as- 
sumed Hubble constant. 

The earliest (optical) studies fou nd that Vtg ry was rel- 
atively constant for 2 < z < 4.5 ([So"ngaila||2001| |Bok-| 
senberg et al.||2003l [Schaye et al.||2003]|!jongaila |2005[ 
Scannapieco et al.||2006p . Cooksey et al.| ( |2010[ ) showed 
that i^civ increased by at least a factor of four from 
3 to z < 1. Using better optical spectra, D'Odorico 
et al. (20101 found that f^civ increased smoothly from 
z = 3 to z = 1.5 and mapped well onto the z < 1 values. 

However, all these studies measured a power-law shape 
for the /(iVciv), which formally corresponds to an infi- 
nite C IV mass density. These surveys were limited to 
small numbers of available sightlines, typically less than 
50. The rarest systems are the strongest, and these dom- 
inate the mass density measurement when the distribu- 
tion is a power law. Hence the small-number statistics on 
the strongest C IV absorbers limited the quoted f^c iv to 
be for 12 < log TVciv < 15, with one sur vey pushing to 
logiVc iv ~ 16 (Scannapiec o et al. 2006). Indeed, scal- 

< 15 



ing the Scannapieco et al.] (|2006|) Ugy to log -/Ve ry 
reduces their value by 45% ( |Cooksey et al.pOlO ) . 
Our analysis, however, provides good statistics o 



on the 

rare, strong systems, which dominate the SDSS sample. 
The observed ficrv can be a pproximated by the s um of 
the detected C IV absorbers (Lanzetta et al. 1991). The 



total column density in a given redshift bin is simply: 

Cl\ (14) 



,2 



a 



Aft, 



2 



E- 2 

A/obs 



where, again, we factor in the counting uncertainty for 
the number of detections crjv obs ■ Then, we estimated the 
mass density relative to the critical density as: 



a 



civ : 



H rn c iVtot 

C(Oc , (Ax(wg) 

2 



(15) 



: ^C IV I 



a N tt 



f <7{AX) 

\{AX(W r )) 



Since the completeness corrections were compiled in 
equivalent width space, we used the median AX(W r ) 
available to each of the absorbers in the given bin, and 
the error was the standard deviation of the detections. 
The small scaling to account for the accepted false- 
positive distribution follows that of the frequency dis- 
tribution (Equation [8]) . 

In Figure |13| we plot fic iv over redshift from several 
studies for togNcr y ^ 14, adjusted fo r the new limits 
and cosmology (see Cooksey et al. 20101 for details) . The 
direct SDSS measurements are lower limits since we have 
predominately saturated absorbers and use the AODM 
to estimate column densities. However, we can still ex- 
clude values below the limits, including the [Songaila 
J2001 1 measurements a t z < 3. 17 This reinforces the 
|D'Udorico et al. (20101 result that ilcw smoothly in- 
creases! In an upcoming paper, we will be combining 



17 |Schaye et ah] ( |2003| and |Songaila| {2005 ) 
depth methods to also measure S Lq jy , but it i 



used pixel optical 
is unclear how best 



to compare their values to our traditional quasar absorption-line 
study. For this reason, we omit discussion of their results. 
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the D'Odorico et al. (2010) and current datasets to fit 
f(f*c iv j and measure TZcrv for log Nc iv > 12 via Equa- 
tion HI 
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Figure 13. Rcdshift evolution of f2civ for logTVrjiv > 14 ab- 
sorbers. The SDSS apparent optical depth values are lower limits 
(blue), since most systems are saturated. However , these lower lim - 
its firmly exclude the z < 3 measurements from |Songaila| ( |2001| ) , 
which have been grayed out with the other superseded measure- 
ments, and support t he result that fir, t v sm oothly increases from 
Z « 3 to z = 1.5 from |D'Odorico et al.| ( |2010} . 



6. SUMMARY 

We have completed the largest C IV surv ey to date 
by leveragin g the SDSS DR7 quasar catalog (Schneider 
et al.||20l6p . In 26,030 sightlines with z QSO > 1.7 and 
(S/N) > 4 pixel" 1 , we identified 14,772 C IV systems 
at least 5000 km s" 1 blueward of the quasar. The entire 
doublet catalog and other tools for analysis (e.g., com- 
pleteness grids) are made available for the community. 18 

The bulk of the C IV survey was automated, from con- 
tinuum fitting to candidate selection. We visually veri- 
fied all doublets in the final catalog. Our Monte Carlo 
completeness tests included the effects of the automated 
algorithms, user bias, and accepted false positives. We 
analyzed the sample as a whole and in ten redshift bins 
of roughly 1500 doublets each. We also constructed and 
analyzed other subsamples as needed, for specific com- 
parisons to published studies. 

The equivalent width frequency distribution f(W T ) is 
described well by an exponential model. The best-fit pa- 
rameters evolved smoothly with redshift, with the largest 
change arising in the normalization. The parameter evo- 
lution follows the trend of increasing dAfcxv I dX with in- 
creasing redshift — a 2.37±0.09 increase over z = 4.55 — ► 
wl.74. For 1.46 <z< 1.96, dhfcw/dX « 0.34. 

In the region of overlap, the SDSS f(W r ) agree well 
with the published column density frequency distribu- 
tions f( Ncw) from smaller, high-resolution high-S/N 
surveys ( |Songaila||2001| JD'Odorico et al.||2010| ). We con- 
verted /(-/Vcrv) to f(W T ) assuming the linear curve-of- 
growth. The published f(Ncrv) distributions were best 
fit by a power-law formalism. Thus, combined with the 



18 See http : //igmabsorbers . info/ 



exponential nature of the strong-end of f(W T ), we at last 
see an exponential cut-off in the equivalent width distri- 
bution that has not been previously characterized be- 
cause of small sample sizes. The location of this cut-off 
sets the high-end range for C IV absorbers and formally 
leads to convergence in the cosmic C IV mass density. 

In this paper, we joined our large sample of strong ab- 
sorbers to the fits of weak absorbers obtained at high res- 
olution. In a forthcoming paper, we will perform a joint 
analysi s of our SDSS sample with the high-resolution 
data of D'Odorico et al. (2010) to fit clean forms to the 
full range ot /(iVcivJ- This will allow us to improve 
our understanding of the C IV mass density evolution, 
where currently our SDSS study only yields lower limits 
for logiVcrv > 14 systems. 

The co-moving line density dVciv /dX smoothly in- 
creases, approximately, ten-fold from z = 6 — > for 
Wr,i548 > 0.6 A absorbers. Physically, this means that 
the physical cross section (<7 p hys) and/or co- moving num- 
ber density of C IV absorbing clouds (n CO m) has in- 
creased steadily over time, since dAfcw /dX is propor- 
tional to their product. We estimated the projected 
C IV-absorbing cross-section of UV-selected galaxies to 
be -Rphys ~ 50kpc, by assuming all the SDSS W T > 0.6 A 
absorbers arise in > 0.5 L* galaxy halos, wh ich is in good 
agreem ent with Adelberger et al. ( 2005 ) and Steidel et al. 

polo] ). 

This is the first paper in a series on various metal lines 
in the SDSS DR7 quasar spectra, which will provide a 
comprehensive and uniform collection of catalogs for fu- 
ture analysis and detailed comparison with cosmological 
hydrodynamic simulations. 
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